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1 Introduction 



Asymptotic heat kernel expansion for nonminimal differential operators on curved 
manifolds in the presence of gauge fields is considered. The complete expressions 
for the fourth coefficient (E4) in the heat kernel expansion for such operators are 
I/-) \ presented for the first time. The expressions were computed for general case of 

manifolds of arbitrary dimension n and also for the most important case n = 4. 
. The calculations have been carried out on PC with the help of a program written 

O ! in C 

o\ 

ON 

a 

B 

>- . In this paper we continue a computer study of the heat kernel expansion of elliptic differ- 
ential operators on curved manifolds and in the presence of arbitrary gauge fields started 
in [0, 0. The coefficients in the asymptotic expansion of the diagonal heat kernel elements, 
called DeWitt-Seeley-Gilkey (DWSG) coefficients [SI EL H], are of fundamental importance 
in quantum field theory, quantum gravity, spectral geometry and topology of manifolds. 
Many quantities of interest (such as the effective action, Green's functions, anomalies in 
quantum field theory 0, ||, |], the indices of elliptic operators and the invariants of 
manifolds in spectral geometry || [|) are expressed in terms of DWSG coefficients ||10|| . 

One more example: the DWSG coefficients, computed for the Sturm-Liouville (1- 



dimensional Schrodinger) operator, form the Korteweg-de Vries hierarchy |TT], [12], |13| . 



As is well known, the heat kernel of a positive elliptic differential operator A of the 
order 2r, acting on a bundle of fc-tensors whose base is a compact closed n-dimensional 
manifold M, has a short-time expansion || || |14| in terms of geometrical invariants 

(x\e- tA \x) ~ Yl E m (x|A)t — , t -> 0+ . (1) 

m>0 



The coefficients E m (x\A) depend on the structure of operator A. Being local covariant 
values of certain dimensionality, these coefficients can be expressed in terms of coefficient 
functions of the operator A, curvature (and torsion, in general case) and their covariant 
derivatives. Up to now the most complete results were obtained for the second-order 
operator 

A = -n + x, (2) 

where □ = g' lu D fl D u , is a covariant derivative including generally different connec- 
tions (affine and spinor connections, gauge fields), X is a matrix in internal space. The 
computation is based on the DeWitt ansatz for heat kernel matrix elements 



oo 

{x\e 

m=0 



*V) = e-^D^^x') Y, E m (x,x'\A)t^, (3) 



where a(x, x') is a geodesic interval which is equal to a half of square of a distance along a 
geodesic curve between the points x and x'\ D(x,x') = g^ 1 ^ 2 det[— D^D v ia\(g')~ l l 2 is the 
van Vleck-Morette determinant, g = detg^. A kernel of the operator exp(— tA) satisfies 
a heat equation. Substitution of the ansatz (f|) into the heat equation leads to recurrent 
relations for the coefficients E m . In the case of the operator (0), the explicit expressions 
for coincidence limits E m (x\A) = E m (x,x')\ x=x > have been computed for m = 0,2,4,6,8 
H, [| [15|, |16| (odd coefficients are equal to zero for closed manifolds, they are reserved, 
however, for manifolds with a boundary). The most complete expression for E$ has 
been obtained by Avramidi [JR3] (a particular case of a scalar field with non-minimal 



conformal connection is considered in |L7|, |18|). In a generalization for the case 
of Riemann-Cartan spaces has been considered. Note that in the flat space case it is 
possible to go much farther in computing the coefficients E m (more exactly trE m ). There 
are corresponding results up to E w pOf []. But even in this case, the limit of computation 
by hands is, apparently, reached already. The natural hope of the further progress is 
connected with development of analytical computer programs for computing the DWSG 
coefficients |T], 0, |23|, |24|, [2"5| . For instance, in pi] the DeWitt algorithm for the operator 



— □ + X in a flat space is implemented on the base of the computer algebra systems 
REDUCE and FORM and in [^] a Mathematica program for computing heat kernel 
coefficients is presented. 

However, the DeWitt method does not apply to higher-order operators and nonmini- 
mal operators whose leading term is not a power of the Laplace operator (for discussion 
see, for example, 0, [^J and a recent paper 0]). The simplest example of nonminimal 



operator is the Navier-Lame operator of classical elasticity, fiAV + (A + /i)V(VVj, involv- 
ing a coupling among the components of a vector-valued function (the Lame constants, A 
and /i, characterize the material). H. Weyl was apparently the first who investigated the 
asymptotics of spectrum of operators of the like type ||29|| . In recent years nonminimal 
operators of a similar sort have been encountered by physicists studying the quantization 
of gauge and gravitational fields in arbitrary gauges |TD[. For example, the quantization 



1 For calculating E\2 (up to total derivative terms) see Refs.gl], §|. We are grateful to C. Schubert 
for bringing our attention to these works. 
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of Yang-Mills field in an arbitrary covariant background gauge leads to the operator 

K> = - 6 ^ nab - - l ) D 7 D ? - 2 f acbG U ( 4 ) 

where □ = D^D^ L is the Laplace operator, is a covariant derivative containing the 
external field potential A^, is the field strength and f abc are the structure constants 
of a corresponding Lie algebra. In an analogous way, the quantization of electromagnetic 
field in an external gravitational field leads to the operator [R{ 



A^ = -g^U - ( - - 1)D„D V + R„ u (5) 

a 

(for an analogous operator in quantum gravity see JJTJ)- Usually operators (f4|), (||) are 
considered in the Feynman gauge a = 1, where they are reduced to a minimal differential 
operator (in particular, ([5]) is reduced in this case to the Hodge-de Rham operator acting 
on 1-forms), and then are treated by the DeWitt method. Some works on computation of 
tiE m are based on the specific dependencies, e.g. of the operator (^), on the parameter 
a [K], [32|]. The nontrivial question worth of study is the interplay between the dependence 
of a heat kernel on the gauge parameter a and various invariants appearing in its expan- 
sion. This in turn is important for calculating the gravitational conformal anomaly for 
gauge fields in a general covariant gauge and for investigating possible gauge parameter 



dependence of the anomaly [32 



The main aim of this paper is to study the heat kernel expansion for a generic operator 
with the structure of @, 

-g^a + aD^D" + X^, (6) 

where □ is now a covariant generalization of the Laplace operator, D 11 is the covariant 
derivative involving affine and bundle connections, X is a tensor field and a is a scalar 
parameter (bundle indices are assumed implicitly). The parameter a must satisfy the 
condition a < 1 in order to ensure the positive defmiteness of the principal symbol and, 
hence, the ellipticity of operator (^). The expansion (|IJ) is valid for nonminimal operators 
as well as minimal ones ||. In view of inapplicability of the DeWitt method we follow 
the approach based on the technique of Widom's covariant pseudodifferential symbolic 
calculus |33[ developed in |34], |35], |36| (for review see also [[37]). At present this is the 



most general method permitting one to handle operators of general type and of arbitrary 
order on curved manifolds. 

There are several computer implementations of the algorithms for second-order mini- 
mal operators based on general purpose computer algebra systems (as is mentioned above, 
in REDUCE and FORM and in Mathematica were used). However, the calcula- 
tion of the DWSG coefficients for nonminimal operators is much more complicated than 
for second-order minimal ones and is out of abilities of these programs for nontrivial orders 
of coefficients. 



The algorithm described in [34] was implemented as a C program [|T], 0]. In 1996 



we rewrote this program. This last version allowed us to get for the first time || the 
complete expression of the coefficient E± for the nonminimal operator (0) (earlier only 



partial results for E 4 were known [10, 38, 39 
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2 General Outline of Algorithm and Its Implemen- 
tation 



In this section we give a short description of the algorithm and the program for computing 
the DWSG coefficients. 

Let us consider a positive elliptic operator which spectrum lies inside a contour C. 
The formula 

allows to express the heat operator exp(— tA) in terms of the resolvent {A — A) -1 of the op- 
erator A. In order to obtain the covariant generalization of the pseudodifferential calculus 
method we use the following representation for the matrix elements of the resolvent 



1 r d n k 

G(x, x', A) = (x\- -\x') = / ^=e u{x ' x '' k) a{x,x',k;X), (8) 

• ! ~ A J (2v^) r^ v /^(a; , ) 



where a(x, x', k; A) is an amplitude, l(x, x' , k) is a (real) phase function which is a biscalar 
with respect to general coordinate transformations, A; is a wave vector. 

The resolvent of an operator A satisfies the equation (A — X)G = 1 which leads to the 
equation for the amplitude: 

(A(x, D a + iDJ) - X)a(x, x', k; A) = I(x, x'), (9) 

where I(x, x') is a transport function having both bundle and Lorentz indices. The co- 
variant generalization of the properties of the flat space phase and transport functions is 
reduced to the following relations [B3 



[{D ai ...D am }l}=0, m>l; [{D ai . . . D a JI] = 0, m > 1, (10) 

where {. . .} means symmetrizing in all indices, and [. . .] means taking the coincidence limit 
(x = x'). Equations Ql0|) together with the "initial conditions" [I] =0, [DJ] = k a and [I] = 
E (E is the unit matrix) allow one to compute the coincidence limits for nonsymmetrized 
covariant derivatives [D ai . . . D am l] and [D ai . . . D am I}. These nonsymmetrized derivatives 
are obtained directly from ( |TUD by reducing all terms to a unified index ordering with the 
help of the Ricci identity. The resulting expressions are universal polynomials in the 
torsion T£ c , curvature tensors Rt cd) W a b and their covariant derivatives. In fact, once 
computed and stored the coincidence limits [D ai . . . D am l] and [D ai . . . D am I] can be used 
in many calculations for different operators A. The functions l(x,x',k) and I(x,x'), 
introduced with the help of formulas (|T0|), play an important role in so-called intrinsic 



symbolic calculus developed in f33[. In fact, just these universal functions manifest the 



geometric properties of a base manifold and a bundle. 

Expanding the amplitude a in degrees of homogeneity of k : a = J2m=i a m(x, x', k; A), 
we obtain the recursion equations for a m from Eq.(^j). For example, for the operator (^J) 
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these recursion expressions take the form 



A ab cr bc 

i ab , 



A ab a lbc + i [-g ab (al + 2D d lD d ) + a(D a D b l + D a lD b + D b lD a )\ a obc = 0, 
A ab a mbc + i [-g ab {al + 2D d lD d ) + a{D a D b l + D a lD b + D b lD a )] <j (m _ 1)6c 
+ (- g ab n + aD a D b + X ab )a {m . 2)bc = 0, m > 2, (11) 

where the matrix 

A ab = g ab( D a lD j _ A ) _ aD a lD b l 

is the principal symbol for the operator (|6]). Solving the recursion equations we obtain 
<j m . The DWSG coefficients are expressed in terms of integrals of the coincidence limits 
\cr m ]: 



E m {x\A) 



d n k 



x [a m ](x, k, A) = J([a m }). 



(12) 



(27r) n y/gJc 2tt 

The integrals in ( |I~2"D can be expressed in terms of the Euler gamma functions and Gauss 
hypergeometric functions for a wide class of operators A. The typical integral takes the 
form 



k ^kai • • • k<i2s 



(k 2r - X) l [(l-a)k 2r -A] f 

9{ ai ...a 2s} ( . j/2n!v( S X n ! 2) jv(l^ ^ +S + n / 2 )/ r ' 1 + m ' ^ 

(47r) n /^2 s rl (ra/2 + s)l (/ + m) 



where g{ ai ...a 2 s} * s a symmetrized sum of products of metric tensors. Using the fact that m 
and I are integer numbers, one can express the hypergeometric function in ([13|) in terms 
of elementary functions with the help of the Gauss relation 

a(l — z)F(a + 1, 6; c; z) = (c — a)F(a — 1, 6; c; z) + (2a — c — az + bz)F(a, 6; c; z), (14) 



and using then the formula [41| 

i)- 1 



\ 1— m 



m-2 



F(l, 6; m; z) 



(1-6) 



m— 1 



m 



m—b—l 



1,2, 



\- (6-m + l) fc fc 
^ jfc! 

m-b^ 1,2,.. 



fc=0 



(15) 



where (a)^ = a(a + 1) . . . (a + k — 1) is the Pochhammer symbol (shifted factorial). During 
simplification of tensor expressions we use various symmetry properties of the tensors 
Rl cd) T^ c) W ab , and also the Ricci identity: 



[D a , D b ]cp 



a\...a k 



/ y L jabra\...a k / < L <nabT a\...ai—\jai+\...a k 



i=l 



i=l 



4- T j D in bl -- bl 4- W un hl - hl 

+ 1 ab 1J L Pa 1 ...a k + ^afe^ai.-.a. 



(16) 
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the Bianchi identities: 



D a R b c de + DdR b cea + D e R h cad + T l ad R b cei + T de R b cai + Tl a R b cdi — 0, 

D a W bc + D b W ca + D c W ab + W ai T* c + W u T m + W c {T\ h = 0, (17) 

and the cyclic identity: 

-"bed + -ttcdb + -"dftc + cd + U c± db + U d^ bc + 1 bi-Lcd + I ci I db + 1 di 1 be ~ U l i& V 

(in this work we don't consider torsion, for computing with torsion see M, Q). 

The above algorithm has been implemented in the C language. The C code of total 
length about 10000 lines contains about 200 functions for different manipulations with 
tensors and scalars. These functions are gathered into two programs COLIM and DWS- 
GCOEF. 

The COLIM program computes coincidence limits of the / and I functions and writes 
them to the disk. Once computed and stored the coincidence limits, being universal 
functions, can be used in many calculations for different operators A. 

The DWSGCOEF program computes E m coefficients by the following steps: 

1. Reading input information (operator, order m, etc.) 

2. Computing a set of asymptotic operators for constructing recursion equations 

3. Computing a m with the help of the recursion equations 

4. Taking the coincidence limit [a m ] 

5. Integrating [o~ m ] to obtain the coefficient E m 

6. Substituting tensor expressions for [D ai . . . D a J] and [D ai . . . D ak I] into E m 

7. Reducing hypergeometric functions to elementary functions in the nonminimal case 

8. Output E m (and its Lorentz trace in the nonminimal case) 

To cut down the intermediate swelling we use term-by-term strategy, i.e., Steps 4-6 
are applied consecutively to single terms of a m generated at Step 3. 



3 Computation of E4 

Note that the coefficient E4 is connected intimately with the dimension of physical space- 
time n = 4 due to the fact that the Atiyah-Singer index of an elliptic operator on a 
manifold of dimension n can be expressed in terms of an integral of E n . We computed 
the coefficient £4 for operator @ on a Pentium-75 PC. For completeness we give first the 
expression for E 2 . In the previous version of the program [0, the result was expressed 
in terms of hypergeometric functions. The current version contains a built-in algorithm 
for reducing hypergeometric functions to elementary functions and is based on formulas 

(S3), (HD- 

E 2 = (4tt)-3 j-CiX" 6 - -j (X ba + g ab X*) - C 3 W ab + C 4 R ab + C 5 g ab R^ , (19) 
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where 



C\ = — 7Tt~ — r-i(l-o) * (4 + An - 6a - 3cm) - 4 - 4n + 6a - 3cm - 2an 2 + an 3 } 
8a(f - 1) 3 1 J 

C 2 = 1 {(1 -a)~% (-4 + 2a + an) + 4 - 2a + an} , 

2a^2 — IJ3 

C 3 = . * rr {(1 -a) 1 "* (8 -an) -8 + 8a-3an}, 
4a(| - 1)2 L J 

C 4 = /w X { (1 - a) - 2 (l2n - 12a - 8an - an 2 + 2a 2 n + a 2 n 2 



24a( 

-Yin + 12a + 8an - 5an 2 } , 

C 5 = - — t~~ — { (1 - aV^ f— 24 + 12a + 8an + an 2 - 2a 2 n - a 2 n 2 ) 

48a(| - 1)3 L v ' 

+24 - 12a - an 2 + an 3 } . 

The full expression for E4 consists of 73 tensor terms with 43 different scalar coefficients^ 
and is given in Appendix. Here we reproduce only the trace (with respect to Lorentz 
indices) tr^E^. The expression for the trace still contains some novelties in comparison 
with the result of Branson et al. who computed tri?4 in arbitrary dimension of a 
space but without gauge field and neglecting the terms with total derivative. 



tr L £ 4 = (4tt)-2 { - C^X* - CzDiDjX 1 ' + C 3 (X/X/ + XyX^ 

•r;.v, ; ,v" + CsXyW* - CeWyX* + rvii',,11'" + r s /.'„/,, 

-CsRijX* - C w R l3 R ij + Cn DR + C 12 R 2 - C 1Z RX* } , (20) 



where 

Ci = ^ -((l-a) 1 ^ (96 - 96a - 48an + 8a 2 n + 6a 2 n 2 + 

-96(1 - a) 2 + 32a 2 n + 2a 2 n 2 - 5a 2 n 3 + a 2 n 4 } 
C 2 = 7T7Z ^((1 -a) 1_ t f_48n + 48a + 24an-4a 2 n + a 2 n 3 



1 


96a 2 (f- 


-2)4 


-96(1 - 


-a) 2 


1 




48a 2 (f- 


-2)4 


+48n - 


48a 


1 




16a(§- 


-1)3 


1 





a 2 n 3 



= 77; /n rr {(1 -a)~§ (-4 + 2a + an) + 4 - 2a + an} 
16a(f-l) 3 ^ J 

C 4 = * > {(l-a)~t(4 + 4n-6a-3an) 

—4 — 4n + 6a — 3an — 2an 2 + an 3 } , 



2 There are 28 linear dependencies among these coefficients. 
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C 5 = 1 {(1 -a)' 1 ^ (l92n + 48n 2 + 576a + 192an- 

\MZCL 

216an 2 - 24an 3 - 1152a 2 - 944a 2 n + 136aV + 62a 2 n 3 + 
2a 2 n 4 + 576a 3 + 664a 3 n + 86aV - 31aV - 5aV - 104a 4 n ) 
-62a 4 n 2 - a 4 n 3 + 2aV - 192n - 48n 2 - 576a - 384an + 72an 2 
+576a 2 + 272a 2 n - 208aV + 10a V + 4a 2 n 4 - 24a 3 n + 26a 3 n 2 
-9a 3 n 3 + aV } , 

C& = loo 2(n_o\ {(1 -a)' 1 ' 1 (-192n-48n 2 + 384a + 816an + 96an 2 

+12an 3 - 768a 2 - 1120a 2 n - 148a V - 8a V - 2aV + 384a 3 
+584a 3 n + 150a 3 n 2 + 7a 3 n 3 - 88a 4 n - 58a 4 n 2 - 5a 4 n 3 + a 4 n 4 ) 
+192n + 48n 2 - 384a - 624an + 48an 2 + 12an 3 + 384a 2 + 304a 2 n 
-164a 2 n 2 + 8aV + 2a 2 n 4 - 24a 3 n + 26a 3 n 2 - 9a 3 n 3 + aV } , 

C 7 = -I -1(1 -a) 1 "? (-96 + 22an + an 2 + 2a 2 n - aV) 

48a(f-l) 2 1 v 1 

+96 — 96a + 26an — 3an 2 + an 3 } , 

(l-a) 2 -f -16 + n 
68 " 180 ' 

C 9 = — {(1 - a)^ (\2n - 12a - 8an - an 2 + 2a 2 n + a 2 n 2 ) 

24a(2—lj3 

-Yin + 12a + 8an - 5an 2 } , 

C10 = — L v- 1(1 -a)~t (_360n + 360a + 296an + 60an 2 + an 3 

1440a(f-l) 3 ^ v 

-112a 2 n - 60a V - 2a V - 4a 3 n + aV ) 

+360n - 360a - 296an + 116an 2 - an 3 + an 4 } , 

C n = -4 1(1 - af-^ ( 480 - 240n - 240a - 120an + 16a 2 n 

480a 2 (f-2) 4 lV ; ^ 

+26aV + llaV + a 2 n 4 + 4a 3 n + 4a 3 n 2 - a 3 n 3 - a 3 n 4 ) 

-480 + 240n + 720a - 360an + 120an 2 - 240a 2 + 104a 2 n - 70a 2 n 2 

+ 15aV-5a 2 n 4 + aV }, 

C12 = } n 7T {(1 -a)-% (-144 + 72a + 56an + 12an 2 

576a(|-l) 3 ^ v 

+an 3 - 16a 2 n - 12aV - 2a 2 n 3 - 4a 3 n + a 3 n 3 ) 
+144 - 72a + 16an - 16an 2 - an 3 + an 4 } , 

Ci 3 = — r-((l-a)~5 f— 24 + 12a + 8an + an 2 - 2a 2 n - a 2 n 2 ) 

48a(f-l) 3 L v ' 
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+24 - 12a - an 2 + an 3 } 



As one can see, the scalar coefficients in (|T9|), (|2~0|) depend in a rather non-trivial way 
not only on the gauge parameter a, but also on the space dimension n. This is a typical 



property of non- minimal and higher order operators | 3~4]| . In the dimension n = 4 the 
coefficients Cj in ( p0| ) take the formQ 

2a -1 . . 8a 2 -21a + 6 ^ a + 4, . 13a 2 + 6a - 24 

Cl = a 2 ln 1 ~ a 5^71 7~> C'a = — ^- ln(l - a) — — — , 

6a 2 36a(l — a) ba z 36a(l — a) 

a 2 _ 13a 2 - 36a + 24 

3 ~ 48(1 -a) 2 ' 4 ~ 48(1 -a) 2 ' 

5a + 8, ,„ . a 5 + 11a 4 + 2a 3 + 74a 2 - 180a + 96 
= — m 1-a ; — , 

5 24a 2 V ; 288(1 -a) 3 

7a- 8, „ v a 5 - 17a 4 + 154a 3 - 362a 2 + 324a - 96 
Cr = —mil — a) ; , 

6 24a 2 V ; 288a(l-a) 3 

1 11 a(4 - 3a) 23a 2 - 46a + 8 

° 7 " 3' ° 8 -~180' ° 9 " 12(1 -a) 2 ' Cl °" 360(1 -a)* 

a 2 - 5a - 2 , , 133a 2 - 168a - 60 

Cii = Kl-a) 



12a 2 v ' 360a(l-a) 
a 2 + 4a - 8 _ 3a 2 - 6a + 4 

12 ~ ~ 144(1 - a) 2 ' 13 ~ 24(1 -a) 2 ' 

In the case of symmetric (with respect to Lorentz indices) matrix X* J the expression for 
Xx^E^ coincides with that obtained in [TO, ESI, [10] (the terms with total derivative were 



not computed there). Note also that in this case the logarithmic dependence on the gauge 
parameter is contained only in the coefficients at the terms with total derivative. 

As to the full expression for E± given in Appendix, we should note its rather cum- 
bersome form. It is quite clear that the next coefficients will be much larger. Thus, 
the problem of conveying computational results of such a kind within scientific society 
becomes rather crucial. In this connection, it seems the proposition to create a special 



electronic archive for results of computations |42| makes sense. Such an archive would be 
useful also for checking by other authors the obtained results and for comparing efficiency 
of different algorithms and programs. For this purpose, however, the problems of iden- 
tification of a canonical basis for tensor invariants, elaboration of efficient algorithms of 
reducing arbitrary tensor expressions to the canonical form implied by such a basis and 
converting from one basis into another one must be solved. 



3 Transition to n = 4 is not a quite trivial procedure because the denominators contain the factors 
n/2 — 2 and this requires to handle uncertainties of the type 0/0. 



9 



4 Conclusion 



The current version of the program computes £4 for the operator @ on a Pentium- 75 PC 
for 4 h 5 min, whereas computation of E 2 (and E 4 for minimal operator) takes trivial time 
(< 1 sec) [|. The main reason for such a tremendous swelling are tensor sub monomials of 
the form 

(A) )ai{A )*• ■ ■ ■ (A Q 1 ) l (A Q ) b . 

For nonminimal operator Aq 1 is the matrix 

/ 1 f Qjkn kh I 

{Aq )ab = -To \" \ 9ab + 



k 2 -\ [ aao (l-a)P- A J ' 

whereas for a minimal operator Aq 1 = fc2r 1 _ A , i.e., a single scalar. However we can cope in 
part with this swelling of intermediate expressions by writing (Aq 1 )^ in terms of projectors 
p3| Piab = Qab — jus and P 2a b = j^-- Using then the standard properties of projectors, 
P™ = Pi, PiPj = 0,i ^ j, we can considerably decrease the computing time. 

Besides, some additional work on structuring large output expressions is needed. In- 
deed, to be understandable or usable, the result of the calculation must be put into some 
kind of standard form. As one can see, any coefficient E m in the expansion (p]) is a linear 
combination of monomials of a given mass dimensionality. These monomials are con- 
structed from tensors by multiplication and index contraction. Factors in the monomials 
could be, in general case, noncommutative and possess symmetries leading to nonobvious 
linear dependencies among the monomials. Thus, we need, first, to identify the linearly 
independent basis of tensorial and invariant functionals of curvature, torsion, gauge fields, 
and other tensors and their derivatives. 

Taking into account additional identities like cyclic, Bianchi etc. makes the problem 
of the choice of a canonical basis rather nontrivial. A certain progress in this direction 



has been reached in j|3], |44[]. For example, in |43[ a partial basis (containing the Riemann 
curvature tensor and its covariant derivatives) for the coefficient E 8 for the operator 
A = — □ + X has been constructed. In |44| an algorithm for constructing a basis for E m 
for the mentioned operator in a flat space has been proposed. The analogous problem 
of construction of a canonical basis for nonminimal operator @ is still waiting for its 
solution. Certainly, the choice of a basis can not be unique and different bases are useful 
for different purposes. 

The next important task is the reduction of an arbitrary tensor expression of a given 
mass dimensionality to a chosen canonical basis for E m . A development of an efficient 
algorithm for such reduction would be an essential contribution on the way to full au- 
tomation of computation of the DWSG coefficients. 

We would like to thank S.A. Fulling from Texas A&M University (College Station, USA) 
for fruitful discussions. This work was supported by Grant No. INTAS-93-2058-ext "East- 
West Network in Constrained Dynamical Systems". 

4 The previous version]!], ^ computed E2 for nonminimal operator for 3 min 15 sec on AT386. 
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Appendix 

In this Appendix we give the full expression for the coefficient E4 in expansion ([!]) for non- 
minimal operator (H) for the cases of manifold of arbitrary dimension n and of dimension 
n = 4. 

E A = (47r)-2 {-C i nX ab -C 2 (nx ba + g ab nXi l + 2D t D a X bl + 2D l D b X la + 2D a D b X l l 

+2g ab D i D j X ij ) - C 3 (^D l D a X tb + DiD b X ai + ^X b W ia - ^WfX**) + C 4 (X^X^ 

+X l a X lb + X a iX bi + X^X* + X b l X la + g ab X l3 X 11 ) + C 5 (X l l X ba + X l a X bl + X b X' la 
+X b X ai + X b t X ai +X ba Xi l + //"''A',' A',' + g ab X i:j X ij ) + C 6 (^'W" 6 + X, a W lb + X b % W ta 

_ w .a x bi _ w b x ia ) + C 7 X a l X lb + c s (x a t w lb - w l a x lb ) - c 9 □ w ab - c w w l a w lb 

-C n Wi b W ia + C l2 W ab XC- Ci 3 DiDjR ijab + CiiDiDjR ta}b — C 15 Ri jk a R ijkb + C ie Ri jk a R ikjb 
+C 17 R ij ab X ij + C I J!,; , '-\V" - C l9 R i a j b X ij + C 2 oRi a j b X ji + C 2l R t ab W 13 - C 22 □ R ab 
+C 23 D t D a R lb + C 2i D l D b R ia - C^R ij R ijah + r, ,,/,',, /,""••''' - C 27 (R t a X lb + R b X ai ) 
-C 28 {R l a X M + R b X ia ) - C 29 (R t a W lb - R b W ia ) + C 3a R t a R zb - C 31 (R^Xf + <i" h l>',, A '■' ) 
+C 32 D a D b R - C 33 RX ab - C 34 {RX ba + g ab RX t l ) - C 35 RW ab + C 36 RR ab + C 37 g ab X tj W lj 
+C 38 //'''U,,A'< + C 39 .'/"''II H '■' + C 40 g ab Ri jk iR m - C 4 i g^RijRV + C 42 g ab □ R 
+C 43 g ab R 2 }. 

The coefficients Cj in arbitrary dimension n: 

Ci = ^ * {(l-a) 1 ^ (-576- 192n + 960a + 432an + 48an 2 - 384a 2 - 288a 2 n - 48a V 

192<r(§-2) 5 L v 

+8a 3 n + 6a 3 n 2 + a 3 n 3 ) + 576 + 192n - 1536a - 336an + 48an 2 + 1344a 2 + 96a 2 n - 96a 2 n 2 
-384a 3 + 152a 3 n + 62a 3 n 2 - 27a 3 n 3 - 2a 3 n 4 + a 3 n 5 } , 



1 



C 2 = T77, r- 1(1 - a) 1 ^ (192 - 192a- 48an + 8a 3 n + 6a 3 n 2 + a 3 n 3 ) - 192 + 384a- 48 

192a 3 (§-2) 5 1 v ' 
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an 



-192a 2 + 96a 2 n - 40a 3 n + 6a 3 n 2 + a 3 n 3 } , 

3 -4 { (1 - a) 1_t (-96 - 48n + 96a + 24an + 48a 2 + 36a 2 n + 6a 2 n 2 - 8a 3 n - 6a 3 n 2 

-a 3 n 3 ) + 96 + 48n - 192a - 24an + 24an 2 + 48a 2 - 84a 2 n - 18a 2 n 2 + 6aV + 48a 3 + 36a 3 n 
-16aV-3a 3 n 3 +a 3 n 4 }, 

* ~r { (1 - a)-^ (-48 - 24n + 80a + 36an + 4an 2 - 24a 2 - 18a 2 n - 3a V) 

64a 2 (§-2) 5 LV v ' 

+48 + 24n - 80a - 12an + 8an 2 + 24a 2 - 10a 2 n - 3a 2 n 2 + a 2 n 3 } , 

^ — { (1 - a)~ 2 (48 - 48a - 12an + 8a 2 + 6a 2 n + a 2 n 2 ) 

64a 2 (f-2) 5 lv v ; 

-48 + 48a - 12an - 8a 2 + 6a 2 n - a 2 n 2 } , 

3 \ {2(1 - a) 1 " % (-48 - 96a- 24an + 72a 2 + 54a 2 n + 9a 2 n 2 - 8a 3 n - 6aV - a 3 n 3 ) 

+96 + 96a + 96an - 336a 2 - 84a 2 n + 18a 2 n 2 + 144a 3 - 4a 3 n - 12a 3 n 2 + a 3 n 3 } , 

-4 { (1 - a)-^ (-48 + 72n + 24n 2 - 48a -156an - 36an 2 + 104a 2 + 78a 2 n + 13aV) 

64a 2 (§-2) 5 lv v ' 

+48 -72n- 24n 2 + 48a + 180an -12an 3 - 104a 2 + 22a 2 n + 53a 2 n 2 - 17a 2 n 3 - 3a 2 n 4 + a 2 n 5 } , 

3~4 { (1 - a) 1-5 (96 + 48n - 96a + 264an + 72an 2 - 432a 2 - 372a 2 n - 90a V 

-6a 2 n 3 + 56a 3 n + 42a 3 n 2 + 7a 3 n 3 ) - 96 - 48n + 192a - 264an - 96an 2 + 336a 2 + 708a 2 n 
+30aV - 36a 2 n 3 - 432a 3 - 228a 3 n + 136a 3 n 2 + 15a 3 n 3 - 7a 3 n 4 } , 

' { (1 - a) 1 ' * (-768 + 1056a + 144an- 240a 2 - 120a 2 n- 8a 3 n- 6a V - a 3 n 3 -4a 4 n 



96a 3 (§-2) 4 

+a 4 n 3 ) +768- 1824a + 240an+ 1296a 2 -456a 2 n + 24aV-240a 3 + 208a 3 n-30aV-a 3 n 3 } , 
1 



{3(1 - a) 1- ^ (128 - 96a + 64an - 128a 2 - 96a 2 n - 16aV + 36a 3 n + 20a : 



n 



96a 3 (§-2) 4 

+a 3 n 3 + 4a 4 n - a 4 n 3 ) - 384 + 672a - 384an + 96a 2 + 720a 2 n - 96a 2 n 2 - 384a 3 - 244a 3 n 
+ 120a 3 n 2 - lla 3 n 3 } , 

1 { (1 - a) 1 "* (-384 + 288a + 96an - 20a 3 n - 12a 3 n 2 - a 3 n 3 - 4a 4 n + a 4 n 3 ) 

4a + 32an + 96a 2 - 48a 2 n + 20a 3 n - a 3 n 3 ) } , 

{(1 - a) 1 - 2 * (96-384a-96an+ 144a 2 + 108a 2 n + 18a 2 n 2 -8a 3 n-6a 3 n 2 -a 3 n 3 ) 



96a 3 (f-2) 4 

+3 (128 - 224a + 32an + 96a 2 - 48a 2 n + 20a 3 n - a 3 n 3 ) } , 
1 

96a 3 (§-2) 5 

+2 (-48 + 240a + 24an - 264a 2 + 6a 2 n + 9a V + 72a 3 - 22a 3 n - 3a 3 n 2 + a 3 n 3 ) } , 

ooon 1 {2(1 - a) 1 -? (-16416 + 58656a + 6072an - 29352a 2 - 16608a 2 n - 966a 2 n : 

2880a J (-2 — 2)4 

+3430a 3 n + 1863a 3 n 2 + 74a 3 n 3 + 296a 4 n - 74a 4 n 3 ) + 32832 - 150144a + 4272an 
+ 176016a 2 - 21504a 2 n - 36a 2 n 2 - 58704a 3 + 16108a 3 n - 738aV - 76a 3 n 3 + 15a 3 n 4 } , 
1 



5760a 2 (f-2) 4 



{ (1 - a) 1 -' (53376 + 35472n - 81840a - 54264an - 6672an 2 + 15992a 2 
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+8766aV + 385a 2 n 3 + 1540a 3 n - 385a 3 n 3 ) - 53376 - 35472n + 135216a + 63048an 
-11064cm 2 - 81840a 2 - 15992a 2 n + 13890aV - 1483a 2 n 3 } , 

55296Qa 4 fn 2 ) { (! - a)- 2 -^ (4534272 + 387072n + 44390400a +20447232an +2022912an 2 

-237126144a 2 - 98031360a 2 n - 10807296aV - 218880aV + 435725568a 3 + 179327872a 3 n 
+21040512a 3 n 2 + 626048a 3 n 3 + 11520a 3 n 4 - 399736320a 4 - 167316736a 4 n - 19284928a 4 n 2 
-44256a V + 61168a V + 712a V + 200829696a 5 + 90282944a 5 n + 9158592aV 
-1199920a 5 n 3 - 237456aV - 7912a 5 n 5 - 58430208a 6 - 32038400a 6 n - 3286544a 6 n 2 
+1093416a 6 n 3 + 257132a 6 n 4 + 14474a 6 n 5 + 10438656a 7 + 7865248a 7 n + 1626552a V 
-189364aV - 103878aV - 9249aV - 625920a 8 - 985792a 8 n - 461792aV - 49696a 8 n 3 
+9512a 8 n 4 + 1508a 8 n 5 ) + 8 (-566784 - 48384n - 6682368a - 2936064an - 277056an 2 
+16842816a 2 + 3514080a 2 n - 564096aV - 105120aV - 14097696a 3 + 886960a 3 n 
+1234272aV - 13600a 3 n 3 - 20376a 3 n 4 + 4928832a 4 - 1609312a 4 n - 180976aV 
+126600aV - 1466aV - 2243aV + 90a 4 n 6 - 1148352a 5 + 239184a 5 n + 44488aV 
-9456aV 3 - 1540aV + 288a 5 n 5 + 45a 5 ra 6 + 78240a 6 + 91000a 6 n + 24600a 6 n 2 - 42040a 6 n 3 
-7155a 6 n 4 + 2160a 6 n 5 + 360a 6 n 6 ) } , 

1 { (1 - a)- 2 - 1 (4534272 + 387072n + 44390400a +20447232an +2022912an 2 

276480(2 ("2 _ 2jg 

-237126144a 2 - 98031360a 2 n - 10807296aV - 218880aV + 435725568a 3 + 179327872a 3 n 
+21040512aV + 626048a 3 n 3 + 11520a 3 n 4 - 399736320a 4 - 167298304a 4 n - 19277248aV 
-48096aV + 59248a V + 520aV 5 + 200829696a 5 + 90209216a 5 n + 9127872aV 
-1184560a 5 n 3 - 229776a 5 n 4 - 7144a 5 n 5 - 58430208a 6 - 31927808a 6 n - 3240464a 6 n 2 
+1070376a 6 n 3 + 245612a 6 n 4 + 13322a 6 n 5 + 10438656a 7 + 7791520a 7 n + 1595832a V 
-1 74004a V - 96198aV - 848 laV - 625920a 8 - 967360a 8 n - 454112aV - 53536a 8 n 3 
+7592a 8 n 4 + 1316aV 5 ) + 8 ( - 566784 - 48384n - 6682368a - 2936064an - 277056an 2 
+16842816a 2 + 3514080a 2 n - 564096a 2 n 2 - 105120aV - 14097696a 3 + 886960a 3 n 
+1234272a 3 n 2 - 13600a 3 n 3 - 20376a 3 n 4 + 4928832a 4 - 1646176a 4 n - 187696a 4 n 2 
+137880aV + 574a 4 n 4 - 2759a V - 1148352a 5 + 239184a 5 n + 44488a 5 n 2 - 9456a 5 n 3 
-1540aV + 288aV + 45a 5 n 6 + 78240a 6 + 91000a 6 n + 24600a 6 n 2 - 42040a 6 n 3 - 7155a 6 n 4 
+2160aV 5 + 360a 6 n 6 )} , 

3 j n {8(1 - a) 1- ^ (24 + 48a - 50a 2 - 16a 2 n + 6aV - 2a 3 - 5a 3 n - 2aV) 

-192 - 192a - 96an + 784a 2 - 16a 2 n - 72aV - 384a 3 + 224a 3 n + 80a 3 n 2 - 28a 3 n 3 
-16a 4 - 36a 4 n + 10a V + 12a 4 n 3 - 3aV} , 

3 1 {(1 - a)- 1 - 2 ? (18432 + 4608n - 55680a - 24576an - 3456an 2 + 64512a 2 

+52992a 2 n + 10176aV + 96a 2 n 3 - 51840a 3 - 48480a 3 n - 11376a 3 n 2 - 528a 3 n 3 + 24a 3 n 4 
+29952a 4 + 21376a 4 n + 5680aV + 488a 4 n 3 -16a 4 n 4 - 5376a 5 - 5680a 5 n - 1756a'V - 92a 5 n 3 
+ 19a 5 n 4 + 512a 6 n +320a 6 n 2 + 16a 6 n 3 - 8a 6 n 4 ) + 12 (-1536 - 384n + 3104a + 896an + 96an 2 
-2272a 2 - 1584a 2 n - 16aV + 88a 2 n 3 + 2048a 3 + 800a 3 n - 440a 3 n 2 - 84a 3 n 3 + 22a 3 n 4 
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+2a 3 n 5 - 448a 4 - 96a 4 n + 28a V + 54a 4 n 3 - 3a 4 n 5 - 192a 5 n + 60a 5 n 3 - 3a 5 n 5 ) } , 

— * . {4(1 - a) 1 "-? (-48n+576a+ 144an - 272a 2 -132a 2 n + 32aV + 12aV 3 -16a 3 
384a 3 (§-2) 5 1 v 

-68a 3 n - 44a 3 n 2 - 7a 3 n 3 ) + 192n - 2304a - 768an + 96an 2 + 3392a 2 - 48a 2 n - 464a 2 n 2 
-24aV - 1024a 3 + 864a 3 n + 408a 3 n 2 - 60a 3 n 3 - 20a 3 n 4 - 64a 4 - 160a 4 n + 4a V + 58a 4 n 3 
-3a 4 n 5 }, 

3 ^ {4(1 - a) 1 -^ (48n + 576a + 144an - 336a 2 - 180a 2 n + 24a 2 n 2 + 12a 2 n 3 - 16a 3 

-68a 3 n - 44a 3 n 2 - 7a 3 n 3 ) - 192n - 2304a - 384an - 96an 2 + 3648a 2 + 144a 2 n - 336aV 
-72a 2 n 3 - 1280a 3 + 800a 3 n + 504a 3 n 2 - 44a 3 n 3 - 28a 3 n 4 - 64a 4 - 160a 4 n + 4aV + 58a 4 n 3 
-3a 4 n 5 }, 

I {(1 - a)- 1 -^ (-8832+3840n+2304n 2 + 13824a-21888an-7296an 2 - 96an 3 

£o\J^±CL y 2" £ ) 5 

+14976a 2 + 30816a 2 n + 9840aV + 768a 2 n 3 - 25344a 3 - 19072a 3 n - 6544aV - 1064a 3 n 3 
-56a V + 5376a 4 + 5680a 4 n + 2332aV + 524a 4 n 3 + 53a 4 n 4 - 128a 5 n - 224a 5 n 2 - 112a 5 n 3 
-16a 5 n 4 ) + 12 (736 - 320n - 192n 2 - 416a + 1872an + 256an 2 - 88an 3 - 1664a 2 - 1088a 2 n 
+360a 2 n 2 + 64a 2 n 3 - 20a 2 n 4 + 448a 3 + 96a 3 n - 28a 3 n 2 - 54a 3 n 3 + 3a 3 n 5 + 192a 4 n - 60a 4 n 3 
+3a 4 n 5 )}, 

3 { (1 - a) 1 - 2 ? (23040 + 11520n + 167424a + 172224an + 32592an 2 - 304320a 2 

-281616a 2 n - 78072aV - 6672a 2 n 3 + 62816a 3 n + 50288aV + 10234a 3 n 3 + 397a 3 n 4 
+6352a 4 n + 1588a 4 n 2 - 1588aV - 397aV) - 23040 - 11520n - 144384a - 172224an 
-38352an 2 + 471744a 2 + 375888a 2 n + 24552a 2 n 2 - 1 1064a V 3 - 304320a 3 - 148496a 3 n 
+35536aV + 8270aV - 1375a 3 n 4 } , 

I {{l-a) 1 -^ (187392 + 41088n - 523776a - 121152an + 2448an 2 + 125760a 2 

+87984a 2 n + 10968a V- 792aV+ 256a 3 n + 376a V + 170a 3 n 3 + 23a 3 n 4 + 368a 4 n + 92a V 



-92a 4 n 3 - 23a 4 n 4 ) + 3 (-62464 - 13696n + 237056a + 22848an - 7664an 2 - 216512a 2 
+33200a 2 n + 12968aV - 1856aV + 41920a 3 - 30160a 3 n+ 1736aV + 1802a 3 n 3 - 263aV) } , 

3 1 {3(1 - a) 1- ^ (-31232 - 8768n + 48256a + 27488an + 3856an 2 - 29632a 2 

( DU(2 Z> J 5 

-27440a 2 n - 7616a V - 652aV + 5536a 3 n + 4496a V + 950a 3 n 3 + 43a 3 n 4 + 688a 4 n 
+ 172aV- 172a 4 n 3 - 43a 4 n 4 ) + 93696 + 26304n - 238464a - 61920an + 1584an 2 + 233664a 2 
+68976a 2 n - 1680aV - 540aV - 88896a 3 - 26096a 3 n + 5152a 3 n 2 + 1154a 3 n 3 - 49a 3 n 4 } , 
1 



{(1 -a)- 2 -% (-61440- 15360a - 34560an+ 1645056a 2 + 379392a : 



2 n 



23040a 4 (f-2) 5 

- 16128a V - 4751616a 3 - 1306752a 3 n - 13632a V + 2400a 3 n 3 + 5413632a 4 + 2036736a 4 n 
+165136aV + 720aV 3 + 476a 4 n 4 - 2663040a 5 - 1520160a 5 n - 258512aV - 17064a'V 
-1384a 5 n 4 + 387648a 6 +458544a 6 n + 154720aV + 22692aV+ 1616a 6 n 4 + 46080a 7 -4384a 7 n 
-28600a V - 9692a 7 n 3 - 884a V - 960a 8 - 7616a 8 n - 3284a 8 n 2 + 644a 8 n 3 + 251a 8 n 4 ) 
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+4 ( 15360 + 34560a + 16320cm - 357504a 2 - 56448a 2 n + 10272a V + 438336a 3 - 4640a 3 n 
-18912aV + 2816a 3 n 3 - 119232a 4 + 38400a 4 n + 2524aV - 4152a 4 n 3 + 305a 4 n 4 + 105aV 
-11040a 5 + 3320a 5 n + 60a 5 n 2 - 170a 5 n 3 + 30a 5 n 4 + 240a 6 - 220a 6 n + 60a 6 n 2 - 5a 6 n 3 ) } , 

Jj- { (1 - a) 1 "* (-27984 - 2472n + 16248a + 9108an + 492an 2 - 1778a 2 n 

1440a 2 (|-2) 4 1 v 

-945a 2 n 2 - 28a 2 n 3 - 112a 3 n + 28a 3 n 3 ) + 27984 + 2472n - 44232a + 2412an + 744an 2 
+ 16248a 2 - 4234a 2 n - 237a 2 n 2 + 91a 2 n 3 } , 

— |(1 - a)- 1 * (192 + 96n - 384a + 72an 2 - 128a 2 - 512a 2 n - 136a 2 n 2 - 4aV 

192a 3 (§-2) 5 1 v 

+416a 3 + 400a 3 n + 118aV + lla 3 n 3 - 56a 4 n - 42a 4 n 2 - 7a 4 n 3 ) - 192 - 96n + 384a - 96an 
-120an 2 + 128a 2 + 656a 2 n + 88aV-44a 2 n 3 -416a 3 -272a 3 n + 146aV + 23a 3 n 3 -9a 3 n 4 } , 

{(1 - a) - * (-192 + 384a - 96an + 64a 2 + 144a 2 n + 32a 2 n 2 - 160a 3 



192a 3 (f-2) 5 

-144a 3 n - 38aV - 3a 3 n 3 + 24a 4 n + 18a V + 3a 4 n 3 ) + 192 - 384a + 192an - 64a 2 
-288a 2 n + 40a 2 n 2 + 160a 3 + 48a 3 n - 34a 3 n 2 + 3a 3 n 3 } , 



1 



., „ - {(1 - a) 1 -^ (192 - 288a + 96an - 144a 2 - 112a 2 n - 20a V + 48a 3 n + 26a' 

tyDQ, y ~2 ~ Z J 4 

+a 3 n 3 + Aa 4 n - a 4 n 3 ) - 192 + 480a - 192an - 144a 2 + 400a 2 n - 52a V - 144a 3 - 144a 3 n 
+70a 3 n 2 - 7a 3 n 3 } , 

{(1 - a) - * (37312+ 12624n+ 1184n 2 - 58976a - 31944an - 51 16an 2 - 204an 3 

960a z (2-— 2) 5 ^ 

+22624a 2 + 21280a 2 n + 6142aV + 599a 2 n 3 + 10a 2 n 4 - 2520a 3 n - 2050a 3 n 2 - 435a 3 n 3 
-20a 3 n 4 - 160a 4 n - 40a V + 40a 4 n 3 + 10a 4 n 4 ) - 37312 - 12624n - 1184n 2 + 58976a 
+ 13288an - 1196an 2 - 388an 3 - 22624a 2 - 1120a 2 n + 2010a 2 n 2 + 85a 2 n 3 - 56a 2 n 4 } , 

2 \ { (1 - a)-^ (-72n + 128a + 96an + 16an 2 - 80a 2 - 68a 2 n - 16a V - a 2 r? 

+8a 3 n + 6a 3 n 2 + a 3 n 3 ) +2 (36n - 64a - 48an + 10an 2 + 40a 2 + 2a 2 n - 7a 2 n 2 + a 2 n 3 ) } , 

* - { (1 - a) 1 "^ (11520 - 2880n + 100416a + 34992an + 2472an 2 - 63072a 2 

5760a 3 (§-2) 5 lv v 

-51240a 2 n - 10836a V - 492aV 3 + 7432a 3 n + 5878a 3 n 2 + 1157a 3 n 3 + 38a 3 n 4 + 608a 4 n 
+ 152a 4 n 2 - 152a 4 n 3 - 38a 4 n 4 ) - 11520 + 2880n - 88896a - 43632an - 1032an 2 + 163488a 2 
+38904a 2 n - 6348aV - 384a 2 n 3 - 63072a 3 - 1072a 3 n + 4862a 3 n 2 - 227a 3 n 3 - 41a 3 n 4 } , 

] n {(1 - a) - * (-288 - 144n + 448a + 224an + 44an 2 + 4an 3 - 112a 2 - 140a 2 n 
192a (-2—2)5 

-56a 2 n 2 - 7a 2 n 3 + 24a 3 n + 18a 3 n 2 + 3a 3 n 3 ) + 288 + 144n - 448a - 80an + 28an 2 - 4an 3 
+112a 2 + 52a 2 n + 16a 2 n 2 - 17a 2 n 3 - 2a 2 n 4 + a 2 n 5 } , 

— * -r { (1 - a)~* (288 - 320a - 96an - 4an 2 + 80a 2 + 68a 2 n + 16a V + a 2 n 3 
192a 2 (|-2) 5 1 v 

-8a 3 n - 6a 3 n 2 - a 3 n 3 ) - 288 + 320a - 48an + 4an 2 - 80a 2 + 20a 2 n - 4a 2 n 2 + a 2 n 3 } , 
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C 35 = 2 ,] l , { (1 - a) 1 ^ (-288 + 144a + 88an + 8an 2 - 24a 2 n - 14a V - a 2 n 3 

-4a 3 n + a 3 n 3 ) + 288 - 432a + 56an - 8an 2 + 144a 2 - 32a 2 n + 14a 2 n 2 - 3a 2 n 3 } , 

C 36 = ^ {(l-a)~% (-432n + 576a + 528an + 144an 2 + 12an 3 - 288a 2 - 392a 2 n 

576a (-^ 2J5 

-176a 2 n 2 - 28a 2 n 3 - a 2 n 4 + 64a 3 n + 64a 3 n 2 + 20a 3 n 3 + 2a 3 n 4 + 16a 4 n + 4a 4 n 2 - 4a 4 n 3 
-a 4 n 4 ) +432n-576a-528an+ 72an 2 -12an 3 + 288a 2 + 104a 2 n + 20a 2 n 2 + 10a 2 n 3 -5aV} , 

C 37 = * {(l-a)" 1 ^ (-192+ 1152a +192an- 1680a 2 - 648a 2 n - 36a 2 n 2 + 768a 3 

+628a 3 n + 96aV + 2a 3 n 3 - 24a 4 - 178a 4 n - 63a V- 5a 4 n 3 - 24a 5 - 2a 5 n + 9a 5 n 2 + 2a 5 n 3 ) 
+192 - 960a - 96an + 720a 2 + 24a 2 n - 36a 2 n 2 - 48a 3 - 20a 3 n + 24a 3 n 2 - 4a 3 n 3 - 24a 4 
+26a 4 n- 9a V + a 4 n 3 } , 

C 38 = I {(l-a)' 1 '^ (-624 - 24n + 1056a + 228an + 12an 2 - 456a 2 - 238a 2 n - 39aV 

192(2 — 2) 5 

-2aV + 24a 3 + 42a 3 n + 21a 3 n 2 + 3a 3 n 3 ) + 624 + 24n - 432a + 108an + 24a 2 - 26a 2 n 
+9a 2 n 2 - a 2 n 3 } , 

C 39 = ^ {(1-a) 1 ^ (576 - 384a- 192an + 44a 2 n + 24a 2 n 2 + a 2 n 3 + 4a 3 n - a 3 n 3 ) 



C40 — 



-576 + 960a - 96an - 384a 2 + 116a 2 n - 4a 2 n 2 - 5a 2 n 3 + a 2 n 4 } 

(1 -a) 2 ~t + n-5 
180(n-4) ' 



C41 = }.„ „, {(l-a)~3 (2160n- 4800a - 3600an - 600an 2 + 3360a 2 + 2984a 2 n + 776a 2 n 2 

5760a^(^-— 2) 5 

+64a 2 n 3 + a 2 n 4 - 448a 3 n - 352a 3 n 2 - 68a 3 n 3 - 2a 3 n 4 - 16a 4 n - 4a 4 n 2 + 4a 4 n 3 + a 4 n 4 ) 
-2160n + 4800a + 3600an - 480an 2 - 3360a 2 - 520a 2 n + 484a 2 n 2 - 54aV - a V + a V } , 



Ci2 = 3 ,» , {(l-a) 1 "^ (l920- 1920a-480an + 64a 3 n + 56aV+14a 3 n 3 +a 

+ 16a 4 n + 4a V - 4a 4 n 3 - a 4 n 4 ) - 1920 + 3840a - 480an - 1920a 2 + 960a 2 n - 320a 3 n 
+64a 3 n 2 - 14a V - a 3 n 4 + a 3 n 5 } , 

C43 = },„ ns {(l-a)-? (1728 - 2112a - 720an - 48an 2 + 672a 2 + 584a 2 n + 152aV 

2304a 2 (|-2) 5 1 v 

+ 16a 2 n 3 + a 2 n 4 - 64a 3 n - 64a 3 n 2 - 20a'V - 2a 3 n 4 - 16a 4 n - 4aV + 4a 4 n 3 + a 4 n 4 ) 
-1728 + 2112a - 144an + 48an 2 - 672a 2 + 104a 2 n - 8a 2 n 2 - 12a 2 n 3 - a 2 n 4 + a 2 n 5 } . 



The coefficients Cj in the dimension n = 4: 

7- 11a + a 2 , , . 84- 174a + 112a 2 - 15a 3 
Cl = 24a3 ln(1 - Q) + 288^(1 -a) < 

-1 + a + a 2 , , -12 + 18a + 8a 2 -7a 3 
° 2 = 24a3 ln(1 - Q) + 288^(1 - a) ' 
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3 + a - 2a 2 , . , 36 - 6a - 36a 2 - 5a 3 
12a 3 Ml~a)+ 144a2(1 _ fl) ■ 

3 , ,„ , 36 - 54a + 12a 2 + 7a 3 
■ ln(l — a) + 



32a 2 v ; 384a(l-a) 2 

1 , N -12 + 18a -4a 2 -a 3 
■ ln(l — a) 



32a 2 v ' 384a(l-a) 2 

1 + 5a - 4a 2 , . . 12 + 54a - 80a 2 + 13a 3 
24a 3 ln(1 a) + 288a 2 (l-a) ' 

13 , , , -156 + 426a- 340a 2 + 71a 3 
■ ln(l - a) + 



32a 2 v ' 384a(l-a) 2 

-3 - 25a + 14a 2 , ,„ . -36 - 282a + 324a 2 + 53a 3 
■ ln(l - a) + 



24a 3 v ; 288a 2 (l-a) 

16 - 18a - 3a 2 + a 3 , , 96 - 156a + 20a 2 + 31a 3 
mfl ln(1 " 0) + 72a 2 (l-a) ' 

-8 - 18a + 30a 2 - 3a 3 , , N -24 - 42a + 121a 2 - 40a 3 
■ ln(l - a) + 



12a 3 v ' 36a 2 (l-a) 

8 - 6a - 6a 2 + a 3 , . , 24 - 30a - 13a 2 + 16a 3 
12^ Ml-a) + 36a2(1 _ a) . 

-l + 7a-2a 2 , , , -12 + 90a - 64a 2 - a 3 
In 1 — a) H — : , 

24a 3 V ' 288a 2 (l-a) 

684 - 2772a + 1863a 2 - 148a 3 , 1368 - 6228a + 6330a 2 - 1409a 3 

360^ ( ~ a) + 720a 2 (l - a) 

-4068 + 4383a - 385a 2 , . , -8136 + 12834a - 3407a 2 
720^ ln(1 - a) + 1440a(l - a) ' 

-1584- 47624a + 31554a 2 + 6483a 3 + 2906a 4 , , 

-. lnfl — a) 

17280a 4 v ' 

+ (-95040 - 2524800a + 11482440a 2 - 18335200a 3 + 12961138a 4 

-3498802a 5 - 762607a 6 + 1727374a 7 - 1316820a 8 + 353710a 9 ) 

1036800a 3 (l - a) 4 ' 

-1584- 47624a + 31554a 2 + 6483a 3 + 3098a 4 , . 

8640^ ln(1 - a) 

+ (-95040 - 2524800a + 11482440a 2 - 18421600a 3 + 13306738a 4 
-4017202a 5 - 417007a 6 + 1640974a 7 - 1316820a 8 + 353710a 9 ) 
518400a 3 (l - a) 4 ' 

-4 - 12a - 9a 2 , , , -24 - 60a - 14a 2 + 83a 3 + 37a 4 
24V ln(1 - a) + 144a 2 (l-a) ' 

-96 + 257a - 101a 2 + 16a 3 



ln(l - a) 



+ - 



288a 3 

-1152 + 5964a - 10458a 2 + 5708a 3 + 1807a 4 - 896a 5 - 1752a 6 + 864a 7 

3456a 2 (l - a) 3 ' 



2 - 10a - 15a 2 , , , 12 - 66a - 62a 2 + 84a 3 + 37a 4 
24^ Ml-a) + I44a 2 (l-a) ' 
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-2-14a-15a 2 , . -12 - 78a - 46a 2 + 102a 3 + 37a 4 
■ ln(l — a) + 



24a 3 v 1 144a 2 (l-a) 

-113+ 173a -40a 2 , , . -1356 + 5466a - 6428a 2 + 1937a 3 - 592a 4 + 1752a 5 - 864a 6 
■ m(l - a) ' 



144a 2 v ; 1728a(l-a) 3 

-180 - 3768a + 4323a 2 - 397a 3 , . . -120 - 2452a + 4158a 2 - 1179a 3 
720^ ln(1 - °) + 480a 2 (l-a) ' 

-916+ 1608a + 39a 2 - 23a 3 , .„ . -5496 + 12396a - 3674a 2 - 559a 3 

■ ln(l - a) + 



720a 3 v ; 4320a 2 (l-a) 

518- 1200a + 1167a 2 - 129a 3 1554 - 4377a + 5042a 2 - 1712a 3 

360^ ( ~ a) + 1080a 2 (l - a) ' 

80 + 520a - 2182a 2 + 1549a 3 - 28a 4 , ,„ 

720^ ln(1 - a) 

1920 + 5760a-87728a 2 +280624a 3 -421348a 4 +331500a 5 -129802a 6 + 18484a 7 +810a 8 + 5a 9 
+ 17280a 3 (l - a) 4 ' 

1578 - 945a + 56a 2 3156 - 3468a + 499a 2 

360^ ( " a) + 720a(l - a) ' 

-3 - 10a + 7a 2 , . . -36 - 66a + 252a 2 - 109a 3 - 38a 4 
m(l — a) + 



24a 3 v ' 288a 2 (l-a) 2 

1 + 2a - 3a 2 , . . 12 + 6a - 68a 2 + 63a 3 - 10a 4 
■ ln(l — a) + 
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24a 3 v ' 288a 2 (l - a) 

-4 - 6a + 13a 2 - a 3 , , , -24 - 24a + 100a 2 - 31a 3 
■ m(l — a) + 



12a 3 v ; 72a 2 (l-a) 

-556 + 355a - 20a 2 , , N -3336 + 7134a - 4327a 2 + 544a 3 
120^ ln(1 - a) + 720a(l - a) 2 ' 

3 - 2a , . , 12 - 26a + 16a 2 - a 3 
W ln(1 " a) + 96a(l - a) 2 ' 

-1458+ 1005a- 76a 2 , , . -2916 + 3468a - 659a 2 
■ ln(l — a) + 



720a 2 v ; 1440a(l-a) 

3 -2a, . , 36 - 46a + 11a 3 
T6^ lD(1 ' fl)+ 192a(l-a) 2 ' 

-3 + 2a, . , -12 + 26a- 16a 2 + 3a 3 
ln(l — a) ' 



48a 2 v 1 192a(l-a) 2 

6 - 7a + a 2 , . . 12 - 20a + 9a 2 
12a 2 h(1 -° )+ 24a(l-a) ' 

9-l2a + 2a 2 , /l , 12 - 34a + 36a 2 - 13a 3 
■ m(l - a) + 



72a 2 v ; 96a(l-a; 
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4 - 28a + 5a 2 , ,„ N 48 - 456a + 988a 2 - 778a 3 + 182a 4 + 13a 5 - a 6 
96^ ln(1 a) + 1152a(l - a) 3 ' 

5a, , , a 2 (180- 450a + 330a 2 -57a 3 + a 4 ) 
— ln(l — a) + 



3 



32 v ' 1152(1 -a) 

-12+ 12a -a 2 , , , -6 + 10a - 3a 2 
24^ 111(1 - a) + 12a(l-a) ' 
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2 - ln(l - a) 
360 ' 

-45 + 60a - 2a 2 -180 + 526a - 452a 2 + 91a 3 

ln(l - a) + 2880a(l - a) 2 

-5 + 5a + 5a 2 - a 3 , , -60 + 90a + 64a 2 - 59a 3 
60^ ln(1 " a) + 720a 2 (l-a) ' 

-9 + 12a -2a 2 -36 + 118a- 116a 2 + 37a 3 

288^ ( ~ a) + 1152a(l-a) 2 ' 
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